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The evolution of the reactor antineutrino spectrum toward equilibrium above the inverse beta-decay threshold 
during the reactor operating period and the decay of residual radiation after reactor shutdown are considered. 
It is found that, under certain conditions, these processes can play a significant role in experiments seeking 
neutrino oscillations. 



1. INTRODUCTION 

The flux and the spectrum of reactor antineu- 
trinos depend not only on the current reactor 
state, which is specified by the power level and 
the isotope composition of the nuclear fuel used, 
but also on the preceding evolution of the fuel. 
The spectrum of antineutrinos produced in the 
beta decay of fission products and other radioac- 
tive nuclei accumulated in the reactor core begin 
to evolve toward equilibrium after the start-up of 
a reactor. After reactor shutdown, the antineu- 
trino radiation diminishes for a long time. 

In [1, 2], equilibration of the spectrum and the 
reduction of the soft section of the v e spectrum, 
where the effects of preceding evolution are pro- 
nounced, were analyzed in connection with the 
problem of searches for the neutrino magnetic mo- 
ment in experiments studying i7 e ,e~ scattering. 

In this study, we consider the inverse-beta- 
decay process 

i7 e + p —> n + e + (1) 

and antineutrino spectrum above the threshold 
for this reaction (1.804 MeV). Presently, reaction 
(1) is of interest mainly as a tool for seeking neu- 
trino neutrino oscillations in reactor experiments. 

In this case, the oscillations in question are 
manifested in the disappearance of some fraction 
of the i7 e flux and in a characteristic modulation 
of the i7 e spectrum, the latter being specified by 
the factor 

P = 1 - sm 2 26sin 2 (1.27Sm 2 R/E), (2) 



where sin 20 is the mixing parameter, 5m 2 is 
the mass parameter measured in electronvolts 
squared, R is the distance (in meters) between 
the source and detector, and E is the detected 
neutrino energy measured in megaelectronvolts. 

In experiments seeking the above oscillations, 
the u e spectrum and flux measured with the aid 
of reaction (1) are compared with their values ex- 
pected in the absence of the oscillations. In this 
case, use is made of the reactor i7 e spectrum ob- 
tained independently (y e spectrum at the produc- 
tion instant). Uncertainties in determining this 
spectrum restrict the sensitivity of the method, 
and systematic error in the spectrum shape may 
generally mimic or mask the oscillation effect. 
Procedures that are applied to analyze data from 
experiments seeking oscillations were described in 
more detail elsewhere (see, for example, [3]). As a 
rule, equilibration of the spectrum and the reduc- 
tion of the v e spectrum do not lead to consider- 
able effects in the region E > 1.8MeV. However, 
their role can be greatly enhanced in some cases 
that will be discussed in Section 4. 

This article is organized as follows. In Section 
2, we give a brief account of available data on the 
spectrum of reactor i7 e in the region E > 1.8MeV 
and on the cross section for reaction (1). In Sec- 
tion 3, we consider the evolution of the spectrum 
in the reactor operating period and residual an- 
tineutrino radiation in the reactor shutdown pe- 
riod and determine the relevant cross sections. 
Section 4 is devoted to discussing the results. 
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2. SPECTRA AND CROSS SECTIONS 
(STANDARD APPROACH 

Let us briefly review the basic properties of the 
reactor-antineutrino spectrum in the energy re- 
gion above 1.8 MeV and the spectrum-averaged 
cross section for reaction (1) (see, for example, [1, 
2, 4-8] and references therein). 

The standard approach is based on the follow- 
ing assumptions: 

(i) The v e spectrum is formed exclusively by 
the beta decays of the fragments produced in a 
reactor upon the fission of 235 [/ ; 2 39 p u 238 jj an( j 

241 Pu isotopes. 

(ii) For each isotope, the equilibrium spectrum 
is established within a small time interval that 
can be neglected. At the instant under consid- 
eration, the reactor-antineutrino spectrum p(E) 
measured in (MeV) -1 units per fission event can 
then be expressed in terms of the spectra p(E) 
for four fissile isotopes. Specifically, we have 

p(E) = J2^iPi(E) (3) 

where on is the contribution of a given isotope to 
the number of fission events occurring in a reactor 
at this instant, the subscript values of i = 5, 9, 8, 
and 1 labeling the quantities associated with the 
fissile isotopes 235 U, 239 Pu, 238 U and 241 F« 

Because of 235 U depletion and the accumula- 
tion of fissile plutonium isotopes, the contribu- 
tions on entering into Eq. (3) change, which 
leads to a change in the total spectrum p(E) over 
the reactor operating period. As was mentioned 
above, the spectra pi{E) are assumed to be time- 
independent. Information about current a, val- 
ues is presented by the reactor personnel, the rel- 
ative error in these values being set to 5%. 

The spectra Pi {E) for 235 U, 239 Pu and 241 Pu 
were obtained in [5] by the conversion method, 
which enables one to reconstruct the relevant an- 
tineutrino spectrum on the basis of the total spec- 
trum of beta- decay electrons from the set of the 
fragments of a given fissile isotope. For a few 
tens of hours, thin layers of the aforementioned 
isotopes were exposed to a thermal-neutron flux 
from the reactor installed at ILL (Grenoble), and 
the current spectra of electrons from the beta de- 
cay of fragments were simultaneously measured 



in the energy region above 2.0 MeV for 235 U and 
in the energy region above 1.8 MeV for 239 Pu 
and 241 Pu . After approximately 12 hours of 
irradiation, the beta spectra reached saturation 
and changed ii longer. Therefore, the antineu- 
trino spectra Pi{E) reconstructed on the basis of 
these (3 spectra are actually those that are estab- 
lished after about one day of fuel irradiation. For 
238 U, use was made of the spectrum calculated in 
[6] because the beta spectrum was not measured 
for the fragments of El isotope. 

The expected number of events of the inter- 
action between antineutrinos and target protons 
was calculated in terms of the cross section ay- a 
for reaction (1), 

ov-A = ^2 ana (4) 

where en — J pi{E)a{E)dE , <r(E) being the reac- 
tion cross section for monoenergetic antineutrinos 
[7]- 

In terms of 10 43 cm 2 units per fission event, 
the Oi values calculated in this way are 

as = 6.39 ± 1.9%, cr 9 = 4.18 ± 2.4% 
ai =5. 76 ±2.1%, cr 8 = 8. 88 ±10% (5) 

These calculations were performed with the9 con- 
stants corresponding to the free-neutron lifetime 
of r = 887.4s ± 0.2%. The resulting error in the 
cross section (4) is 2.7% (68% C.L.). 

To a higher precision, the cross section for reac- 
tion (1) is known from an experiment performed 
by the Kurchatov Institute - College de France - 
LAPP collaboration at a distance of 15 m from 
the Bugey PWR reactor [8]: 

Cmeas = 5.750 x 10~ 43 cm 2 / (fissionevent) 

±1.4%(68%C.L.) (6) 

This result corresponds to the following contribu- 
tions a; from fissile isotopes: 

a 5 = 0.538, a 9 = 0.328, a 8 = 0.078, 

ai = 0.056. (7) 

The proportion in (7) for the number of 
isotope-fission events is typical for PWR reactors, 
which were used in the majority of experiments 
seeking the oscillations in question. Within the 
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Table 1 



Yields of some fission fragments (in %). 



Fission 
fragment 


Multicolumn4Fissile nucleus 




235 jj 


239p u 


241 Pu 


238 




5.95 


5.30 


4.89 


5.50 


132j 


4.30 


5.40 


4.14 


5.16 


93 y 


6.40 


3.89 


3.15 


4.97 


106 Ru 


0.40 


4.31 


6.18 


2.55 


144 Ce 


5.48 


3.74 


4.39 


4.50 


90 Sr 


5.82 


2.10 


1.57 


3.12 



errors, the cross section in (6) agrees with the 
reaction cross section oy-A found for the given 
composition of nuclear fuel: 

Cmeas/cv-A = 0.987 ± 1.4% (experiment) 

±2.7%(V-A). (8) 

The experimental value in (6) is treated as a 
metrological reference for the cross section in the 
absence of oscillations. A feature peculiar to this 
reference is that, in each specific case of its appli- 
cation, it must be rescaled to the relevant com- 
position of nuclear fuel. As a result, the error in 
the cross section increases to about 1.6%. 

For reactors of the type being considered, the 
duration of the operating period is approximately 
one year. After that, the reactor is shut down for 
30-40 days, and one-third of the fuel is replaced by 
a new load of fuel. Therefore, the fuel is irradiated 
for three years. 

The detector background is measured during 
reactor shutdown periods. It is assumed that an- 
tineutrino emission in the region E > 1.8 McV 
ceases completely within one day after reaction 
shutdown. 

The above information forms a basis for ana- 
lyzing and interpreting the results of experiments 
seeking neutrino oscillations in reactor experi- 
ments. 



3. SPECTRA AND CROSS SECTIONS 
(INCLUSION OF NONEQUILIB- 
RIUM EFFECTS) 

1. Let us determine more precisely sources 
that contribute to the formation of the rcactor- 
antineutrino spectrum p(E), 

p(E) = F p(E) + c p(E) (9) 

The first term in (9) describes radia- 
tion from the set of fragments produced in 
235 [/, 239 Pu 238 U and 2A1 Pu fission, their inter- 
actions with neutrons being disregarded here; the 
second term takes into account additional radia- 
tion arising in radiative neutron capture by accu- 
mulated fragments. 

For an irradiation time of t on > 1 day, each 
of the spectra p(E, t on ) for the four isotopes 
235 U 239 Pu 23S U and 2A1 Pu were calculated as 
functions of the duration of fuel irradiation. In 
this calculation, the fission rate was assumed 
to be constant for each isotope. The residual- 
radiation spectra p(E,t ff) were also calculated 
as functions of the time t a f / from the end of fis- 
sion. 

The calculations were performed by summing, 
at each instant, the contributions from individual 
fission fragments with allowance for their yields 
from the fission process, decay diagrams, and life- 
times. The database that we used contains infor- 
mation about 571 fission fragments whose cumu- 
lative yields are not less than 10~ 4 % each. As a 
matter of fact, the activity of the majority of the 
fragments that can contribute in the energy re- 
gion E > 1.8 MeV reaches saturation within one 
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day after the onset of the fission process, and a 
further increase is due to only six fission products 
whose properties are well known. Three of them 
- 97 Zr {E max = 1.922 MeV) 132 I {E max = 2.140 
MeV), and 93 V (E max = 2.890 MeV) - attain 
equilibrium within ten days. A further slow in- 
crease is determined by the 106 Ru and 144 Ce half- 
lives (see (10)). 

Finally, some contribution comes from 
9 "Y(T 1/2 = 64 h, E max = 2.279 MeV) as well, 
which is in equilibrium with its very long-lived 
predecessor 90 Sr (Ti/ 2 = 28.6 yr). The yields of 
these fission fragments are quoted in the table. 

2. The antineutrino spectra calculated for 235 U 
and 239 Pu are displayed in Figs, la and lb (equi- 
libration) and in Figs. 2a and 2b (decrease). For 
four fissile isotopes, Figs. 3 and 4 show, respec- 
tively, the cross sections <jd n) for reaction (1) 
in the fission process and the decrease in <J(t ff) 
in the residual-radiation spectra. In Figs. 1-4, 
the spectra and cross sections are presented in di- 
mcnsionless units and are normalized to the cor- 
responding values after a lapse of t on — 1 day 
from the onset of the fission process. 

3. The correction p(E) was calculated with al- 
lowance for fragments accumulated in a reactor, 
the cross sections for radiative neutron capture, 
and the spatial and energy distributions of the 
neutron flux over the reactor core. The results 
presented in Fig. 5 refer to a standard operating 
period of a PWR reactor. 

The relevant contribution ° p{E) to the cross 
section for reaction (1) is about 0.2% at the end 
of the operating period. After reactor shutdown, 
the spectrum c ' p{E) decreases fast, not making a 
significant contribution to residual radiation. 4. 

4. DISCUSSION OF RESULTS 

(i) In the approximation specified in Section 2, 
the flux and the spectrum of reactor antineutri- 
nos are un- ambiguously determined by the cur- 
rent reactor state. After sharp changes in this 
state, the characteristics of the v e flux take the 
corresponding equilibrium values within a time 
as short as one day. In particular, the antineu- 
trino flux falls down to zero within one day after 
reactor shutdown. 



As a matter of fact, the antineutrino flux in 
the region E > 1.8 MeV has been found to have 
a noncquilibrium component, whose relaxation 
time exceeds the duration of the reactor opera- 
tion period. In a view of this, a determaination 
of the i7 e spectrum and of the corresponding val- 
ues of the cross section for reaction (1) requires 
tracing the preceding evolution of reactor oper- 
ation over a long time and taking into account 
power levels, shutdown periods and the discharge 
of spent reactor fuel. In each specific case, this 
can be done, if needed to a sufficient precision. 

The sample of the results in Section 3, which 
were obtained under assumption that the fission 
of uranium and plutonium isotopes proceeds at a 
constant rate, has enabled us to reveal qualitative 
features of effects induced by the nonequilibrium 
component. 

(ii) First of all, wc note that the resulting cor- 
rections to the spectrum and cross sections are 
relatively small, but they are not negligible. In 
the antineutrino energy range 1.8-3.5 MeV, the 
relative contribution of the additional radiation 
during the reactor operating period (Figs. 1, 5) 
is about 4%, which is somewhat greater than the 
error of the ILL spectra [5] . 

The corrections to the cross sections cr, in Fig. 
3 and the correction associated with radiative 
neutron capture may change the cross sections 
<jv-a by 0.4-0.6%. As a result, relation (8) will 
change accordingly. Corrections on the same or- 
der of magnitude may arise if the cross section (6) 
is used as a reference value for the cross section 
in the absence of the oscillations. 

(iii) Here, we consider a situation where resid- 
ual radiation from a stopped reactor can play a 
significant role and provide numerical examples 
illustrating the scale of the effects under discus- 
sion. 

Let us consider an experiment where antineu- 
trinos from reactors are recorded by one detector 
positioned in such a way that the reactors are at 
markedly different distances from the detector. 
Such an experimental setting was implemented, 
for example, in Rovno (distances of 18 and 98 m) 
and in the Bugey-3 experiment (15 and 95 m) 
[9] ; three reactors in Krasnoyarsk were located at 
distances of 57, 57, and 234 m from the detector 
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ma Ti in = 372day Ti / 2 = 30s 

W6RU E±^0MMeV - Rh E max UuiMeV " ^ 

E max = 0.32MeV E max = 2.996MeV v 1/2 7y v ' 



[10]. The idea of these experiments consists in 
using the shutdown and operating periods to de- 
termine the background level and the signals from 
each individual reactor. By comparing these sig- 
nals, one can reveal the oscillation effect or set 
limits on the oscillation parameters. 

By way of example, we consider two identical 
PWR reactors of thermal power 2.8 GW each and 
a detector positioned at distances of 15 and 100 
m from the near and the far reactor, respectively, 
so 

Fig. 5. Ratio of the spectrum component ?(A) 
associated with neutron capture by fission prod- 
ucts in a PWR reactor to the fission spectrum 
p(E). Curves /, 2, and 8 correspond to the be- 
ginning, middle, and end of the reactor operating 
period. 

that the signal from the near reactor is approx- 
imately 45 times as great as that from the far 
one. In this case, the numbers of events of re- 
action (1) that are induced by the near and the 
far reactor per 1 t of a Cffi.s target per day are 
N°P ear » 12000 day- 1 t" 1 and N° p ar » 260 day— 
t _1 , respectively. 

When both reactors are shut down, the detec- 
tor records, however, the residual interaction in 
addition to the background. Within one to two 
days after the shutdown of the near reactor, the 
number of interaction events induced by the resid- 
ual radiation from it in the detector per unit time 
decreases by a factor of about 200, falling down 
to a value of N°H r « 70 day— t~\ which is 
approximately 25% of the signal from the far re- 
actor. After that, the signal from the residual 
radiation decreases smoothly (see Fig. 4) until 
the discharge of the spent fuel begins, which sig- 
nificantly affects detector readings. 

Obviously, these cases require a new approach 
to processing and analyzing experimental data. 



5. CONCLUSION 

The effect of the equilibration of the reactor- 
antineutrino spectrum on the cross section for the 
inverse-beta-decay process (1) with the threshold 
of 1.8 MeV has been considered. The increase 
in the antineutrino flux due to an increase in the 
cross section (this effect was previously ignored) 
is about 0.6% over the reactor operating period. 
This value is commensurate with the accuracy 
(the standard error is 1.4%) in the measurement 
of the cross section itself. A correction at this 
level must also be introduced in the ratio of the 
measured cross section to the cross section calcu- 
lated on the basis of the antineutrino spectrum 
that was determined independently and which 
corresponds to one day of reactor operation. 

We have calculated the additional contribution 
to the antineutrino spectrum from neutron cap- 
ture by fission products and determined the cor- 
responding increase (0.2%) in the cross section for 
reaction (1). It has been found that corrections 
previously ignored in the standard experimental 
scheme aimed at searches for neutrino oscillations 
(one detector and near and far reactors) may be 
as large as 25%. 
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Figure 1. Ratio of the antineutrino spectra for 
(a) 235 U and (b) 239 Pu fission to the spectrum 
after one day of irradiation. The numbers on the 
curves indicate the irradiation time. 
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Figure 3. Ratio of the inverse-beta-decay cross 
sections for fissile uranium and plutonium iso- 
topes irradiated for the time t on to that after one 
day of irradiation. 



Figure 2. Ratio of the antineutrino spectra of the 
rcsidiul radiation from (a) 235 U and (b) 239 Pu ir- 
radiated for two years to the spectrum after one 
day of irradiation. The numbers on the curves 
indicate the time that have elapsed after the ter- 
miation of irradiation. 
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Figure 4. Ratio of the cross sections for inverse 
beta decay induced by the residual-radiation an- 
tincutrinos for uranium and plutonium isotopes 
irradiated for two years to the cross section after 
one day of irradiation as a function of the time 
t a f f after the termination of irradiation. 
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Figure 5. Ratio of the cross sections for inverse 
beta decay induced by the residual-radiation an- 
tincutrinos for uranium and plutonium isotopes 
irradiated for two years to the cross section after 
one day of irradiation as a function of the time 
t ff after the termination of irradiation. 



